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ABSTEACT 

The Bake troposcatter ■bechnique as a tool for cha- 
nnel sonnding has been studied in this report with reference 
to the Kanpur - Hainital tropo scatter link, The technique 
consists of transiaitting a pseudo randon (PH) sequence of 
pulses over a troposcatter channel. At the receiver, the 
received signal is correlated with identical replicas of, 
the PN sequence, generated locally, with different delays 
introduced to account for the multipath delays. Such a 
technique yields delay - doppler plots tezmied scattering 
functions which characterisa the channel. 

Digital computer programs are developed through 
this, study to simulate the results of Eake troposcatter 
technique based on a certain model for scattering. The out- 
put of a single correlator, called a tap yields a doppler 
spectrum which depends upon a number of factors. Scattering 
volume available with in the mid path plane region common to 
a scattering layer and a choosen time delay shell serves 
as a first estimate of the shape and location of the doppler 
spectrum. Antenna gain patterns, antenna off sets, angular, 
dependence functions and Bake modulation signal ambiguity 
functions are employed as modifiers or weighting functions 



of the first estimate to obtain, a final estimate of the 
doppler spectrum for the choosen delay. 

Ibe scattering functions for various beam swinging 

experiments, for different, scattering layer heights, sca~ 

ttering layer thickness and for multiple layers having 

different degree of anisotropy are simulated. The relative 

zero doppler power received at each Bake tap, which gives 

the relative measure of retractive index structure function 
2 

coefficient 0^ , is evaluated. 

The resulting scattering functions indicate the role 
of various factors such as heights of scattering layers 
and their thicknesses, wind velocity, anisotropy of sea-* 
ttering turbulence, wave-number dependence of the spectrum 
of turbulence, etc,, on the receiver output. It is shown 
that these parameters can be studied with the help of a 
Bake system. 
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INTROD-trOTIOH 

Background 

1.1 Brief Hi.storlcal/of Hake Concent and MeasureEents, of 
Bcatterinsr Btinction 

Tlie Bake concept of neasurenent involves transmi- 
ssion of a special wide - bandwidtli signal through, the 
tropospheric scatter channel to a correlation type of radio 
receiver. The signal, in its passage through the scatter 
cliannelj, encounters a dynamic multipath situation- that 
rather uniquely disperses (or spreads) the signal in both 
tine delo-y and frequency. The receiver output signals can 
be employed through further computer processing for obtai- 
ning delay doppler power spectruns or scattering functions. 
The scattering function graphically display the doppler 
spectrum of the signal received with different dolaj/’s.. 

The Bake concept was first introduced by Price 
and Green [l] in 1958 for the purpose of solving the mul- 
tipath problems associated with high frequency digital 
data communications via the ionosphere. Only in more re- 
cent years has the concept been applied to measurement of 
tropospheric scatter propagation characteristics by 
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Barrovr [2] and Birkoneier E3]. Barnow [2] carried out 
the experinents over two paths - a 480 kns overland path 
and a 614 kns over water path. The ain of his experi— 
nents was to record the randonly varying in phase and 
quadrature conponents of each Bake tap output and thus 
investigate the fluctuations in signal phase, as well 
as in anplitude caused by fading. Scattering functions 
were obtained, which clearly illustrate the distribution 
of received signal energy as a function of nultipath tine 
delay and doppler frequency shift, Birkeneier [3] carried 
out experinents for remote sensing of the lower atnos— 
pheric refractive structure. The detected atmospheric 
structures showed spatial refractive inhonogeneity in 
the form of discrete scattering layers at different hei- 
ghts, Birkeneier [4] in 1975 carried out Hake soundings 
over a pathlength of 270 kns and the results have been 
compared with mat eoro logical soundings of the refractive 
index structure determined from radiosonde "ascents 
near the common volume. It was observed by him that the 
received signal zero - doppler power profile indicated 
more power than was predicted by refractive index stru- 
cture function coefficient , 
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1.2 Details of the Problen 

A tropo scatter link has been set-up between 
Eanpur and Uainital and a Rake sys ten hardware is being, 
developed for this link. The data obtained from the 
Rake systen can be usefully enployed to study sone tro- 
pospheric characteristics. In this report sone theoriti- 
cal calculations are carried out to indicate the type of 
output expected fron the Rake systen under various neteoro- 
logical conditions. 

A conputer progran has been developed to plot 
scattering functions, for the paraneters of the experi- 
nental Kanpur - Rainital troposcatter link using a Rake 
systen. Scattering functions are sinulated for the fo- 
llowing connon volune characteristics and experinental 
set-ups i ■ 

(a) Single scattering layer in the connon volune 

(b) Multiple scattering layers in the connon 
volune 

(c) Antenna azinuthal and elevation off sets 

(d) Different values of power index ’p' of 
power spectrun of refractivity fluctuations. 

The relative zero doppler power received at 
each Rake tap, which gives the relative neasure of 
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2 

Refractive Index Structure Function coefficient On , 
is evaluated, 

1.5 Outline of Work 

This work has been divided into various Chapters. 

In Chapter 2;^ link details, connon volume geonetiy and the 
modulation signal to be used have been given in brief. 

In Chapter 5» basic Rake neasurenent concept is given in 
detail, so that one can visualize how the concept is app- 
lied in the development of simulation program. Experi- 
mental Rake System for Kanpur - Nainital troposcatter 
link is also discussed in this chapter. In Cha.pter 4 
tropospheric scattering model as described by Birkeneier[3] 
is given in conjunction with other areas that effect the 
received signal characteristics such as antenna off set 
effects, angular dependence, ambiguity function and 
magnitude weighting function. In Chapter 5 n digital 
computer simulation program is developed based upon the 
Rake concept of measurements as developed (in Chapter 3) 
and the tropospheric scattering model as given (in 
Chapter 4), Then scattering functions for various co- 
mmon volume characteristics and the beam swinging 
experiments are simulated and plotted. Chapter 6 treats 
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2 

the refractive index stiructure function coefficient Cn 
and its evaluation by the zero doppler power profiles. 
Chapter 7 deals with several applications of the digital 
conputer siaulation of scattering functions and con““ 
eludes this report. 

Siaulation progran has been written in Fortran 17 
language and IBM-1800 is used for this work, Scattoring 
Functions and constant delay shells have been drawn by 
IBM-1627 plotter. 
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CHAPTBE 2 

EA-EPUn - HAINITAI TSOPO SCATTER LINK DETAILS 

2.1 linls: Paraneters 

The Kanpur - ETainital troposcatter link has been 
described by Gupta [5] and its geonetry is shown in 
Fig, 2,1, The frequency of operation of the link is 

2.1 GHz (a wavelength of 14.286 cns). Transmitter and 
receiver antennas are situated at a height of 4.3 Mtrs 

and 18,4 Mtrs above ground level (1892,3 Mtrs and 145.8 Mtrs 
AMSI respectively) at Hainital and Kanpur respectively. 

The effective radius of earth is taken as 8500 KMS, Total 
distance (2D) between transmitter and receiver is 330 KMS. 

2.2 Common Volume Geometry 

Calculations in respect of common volume geometry 
have been done by Gupta [5] and the same are used for 
simulation. The following are the common volume parame- 
ters which are used in this work : 

(a) Grazing ray height (Z^) = 1,8 KMS, 

(b) Scattering angle (©) = 38,8 MR, 

(c) Length of the common = 141 KMS, 


volume 




P'G.M link geometry 





FIG. 2.2 PLAN AND ELEVATfON VtEWS OPf A 
TROPOSCATTER LINK [THE LIMttSfW 
AND WIDTH OF SCATTERING VOLUME A|E: ' 

Taken at 3db potNTsI 
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(d) Height of the centre of the connon 

volume above earth’s surface =3.2 KMS 

(e) Diagonal width of the comon 

volume = 2,8 EMS, 

Comon volume geometry for this tropo link is shown in 
Dig. 2,2. 

2 , 3 Antenna Gain Pattern 

The product of transmitter and receiver antenna 
gain is referred to as combined antenna gain and can be 
expressed as 

G = (X,Y,2) Gj^ (X,T,Z) (2.1) 

Throughout the simulation work in this report Gaussian 
circular antenna gain pattern is assumed. The combined 
antenna gain in T~direction by using Gaussians can be 
expressed as 

G = Exp [-. (Y - Y^)^/Y^^] (2.2) 

where, 

Y = Cross path distance in EU'K from the centre of 
the beam 

T = Bean off set distance in KMS at the nid path 
c 

P 2 

Y = Variance (2.6 EMS for the antenna used on 
m 

Kanpur Hainital troposcatter link). 
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The radius of the circle of the antenna pattern at 
the nid path is calculated as 3 hns for a 3 db bean width 
of 0,96 degrees. 

2,4- Modulation SiCTial 

"Very short periodic pulses spaced sufficiently 
apart so that the nultipath response die out between succ- 
essive pulses, can be used for channel sounding. The 
disadvantage is high peak to average power ratio. In a 
Bake systen this is overcone by using a constant anpli- 
tude signal - a radio frequency carrier that is pseudo- 
randonly nodulated in phase by a PN sequence obtained fron 
a PN sequence generator. PK sequences are choosen because 
(i) their auto correlation function has a central peak 
of value equal to length of sequence and the side-lobe 
is of unit height, (ii) the cross correlation of a PM 
sequence and its shifted version is zero. 

On the Kanpur - Mainital troposcatter link, at the 
transnitter, a PK sequence of length (2'^ - l) i.e., 

32, 767 bits is used to phase modulate the carrier at 
10 M bauds. So the transnitted, PK sequence period (T) 
is 3.2767 nilli secs. At the receiver, the received 
signal is split into in phase and quadrature components 
which are cross correlated with identical binary streams 
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derived fron PIT sequence generators, each of which is 
delayed successively by 0.1 nicro secs. This provides 
a nultipath delay resolution (t) of 0,1 nicro secs^ 
However, a nultipath delay resolution of 0.05 nicro-secs 
is used for sinulati-on in this work. 
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CHAPTER 3 

BRIEF BISCUSSIOIT ON RARE SYSTEM MEASUREMENTS 

3.1 Basic Bake Measurenent Concept 

Rake as a comimication systen is designed to 
work against the combination of random multipath and 
additive noise disturbances. This has been proposed by- 
Price . and G-reen [l]. In channel sounding experiments, 
using Rake principle, the interest is to get a continuous 
record of quadrature components of channel imp-ulse res-’ 
ponse p(T,t) as suggested by POOL [6], From this data, 
it is possible to observe multipath structure, doppler 
shifts and calculate the scattering function which ade- 
quately characterise the tropo-channel. 

As an over-the-horizon propagation medi-um, the 
troposphere is ass-umed to have the form of a linear time 
varying filter. So the input - output relationships 
of the tropospheric channel may be specified by the 
equivalent low pass impulse response p('c,t) which is 
output response of the channel, measured at time t, due 
to a unit impulse input applied at tine (tj-sc)*. The term 
* Rake • applies to a class of measurement methods wherein 
a transmitted signe-d is modulated in a manner- that allows 
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the received signal nodulation to "be cross-correlated 
with delayed replicas of the transnitted signal modu- 
lation to provide estimates of ^(T,t) at several speci- 
fic tine delay values, where k = 1, 2, — : — N. 

Let uo(t), be a series of unit impulses eq.ually 
spaced in time at intervals of T, starting at a time to, 

i.e* uo(t) = ^ — 6[t— (to + nl)] (3*1) 

n 

n = 0,1,2, 

If u(t) be the received signal, then u(t) is given by 
a.> 

u(t) = ^ P(T,t) uo(t - t) dT (3«2) 

-cc- 

which is the convolution of the channel input with the 

impulse response of the channel* Substituting (3*1) i^^ 
(3*2) and interchanging the order of summation and inte- 
gration, u(t) is given by 

f I 

~ <' J “ (to+nT)— di; ^ (3*3) 

n j_-co j 

Using the shifting properly of the unit impulses 

u(t) = 4-nT);t] (3*4) 

n 
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Fron (3*4) it is observed tkat the received signal is 
a series of impulse responses resulting from unit inpul 
ses applied at time (to + nf) . 

To calculate scattering functions, it is essen- 
tial to measure p(t, t) at fixed values of delay, i,e* 
measurements of t). t) can be obtained by 

cross-correlating the received signal with a delayed 
replicas of the transmitted signal, as given by 

. r 

Ruuo = Y J ^(■^) iio^(t - Tg.) dt (3.5) 

I 

Substituting (3*1) and (3.4) in (3.5)* 

Ruuo (t^) = Y p[t-(to + mT) ,*t] 

6[t - (to + nT) -rg] dt (3*6) 

n 

The value of (3.6) is z-ero, except when n = n. Inter- 
changing the order of summation, and integration, 

Ruuo (vg.) = Y p(t-(to+nT) f t] 

n I 

6[t - (to + nT) - Tj-] dt (3*7) 

where I = integration interval =1. 
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RuudCt^.) = “ = I (3.8) 

The neasurenent process involved in obteirmig the equation 
(5.8) forn the basis of the Rake concept which is illu- 
strated in Fig. 3.1. 






Figure 3.1 

The structure of pCir , t) can be estinated by paralleling 
a sufficient number of tap units, each having a different 
reference delay. Consecutive tap unit outputs form a 
series of unit impulses representing the auto correla- 
tion of the transmitted modulation signal at specific 
values of T, A plot of such a series of unit impulse 
auto— correlations resembles the teeth of garden rake 
in appearance, providing a basis for the term ’Rake'. 

3.2 Tap Unit Spacing Considerations 

Minimum tap unit spacing is a function of the 
time resolution capabilities of the transmitted modulation. 
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or probing signal, Maximun tine, or delay resolution 
obtainable with a given probing signal is an inverse 
function of the effective bandwidth of the signal, 

1 

i . e , , Tap unit spacing = ■ 

Effective bandwidth 

If , Tap unit spacing = t , and effective bandwidth is 
WHz, then. 



For the experimental set up of Kanpur - Hainital 
troposcatter link a Rake system is being built with mul- 
tipath resolution of 0,1 micro -secs. For the simulation 
work, however, a delay resolution of 0.05 micro secs is 
used. There are. 19 Rake taps, delays from 0,05 IJ-sec -to 0,95 psev 
with increments of 0,05 psec. 

Rake System for KanoTur - Nainital Troposcatter Link 

The transmitter in a Rake system consists of the 
following sub -systems : 

(a) a highly stable oscillator from which other i 

desired freq_uencies are synthesised 

(b) PR sequence- generator 
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(c) phase modulator f and 

(d) band pass filter 

The receiver consists of s 

(a) highly stable oscillator (synchronised with 

the transmitter oscillator) and frequency 
synthesiser ,, 

(b) Mixer 

(c) IF Amplifier 

(d) Demodulator 

(e) PN sequence generator 

(f) Two correlators at each Rake tap one for 

estimating inphase and the other for quadra- 
ture component of t). 

At the receiver the received signal is demodulated to 
get the inphase and quadrature components. This can be 
achieved by using the arrangement given in Fig, 3*2. 



Inphase 

component 


Quadrature 

component 


Fig* 3.2 
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For the experimental set up on Kanpur - Nainital 

tropo scatter link two Bake systems are considered based 

on the Bake systems built at University of Wisconsia 

(Rake I_jBirkemeier 1969) and at S^lvania (Rake II, 

Rake I and 

Barrow 1965). The transmitter is sane for both^/Rake II 
as shown in Fig. 5.3. 

Rake ^stem I employs analog signal processing with 
attendant requirements on gain and frequency stabilities 
and is expected to have large dynamic range. Rake Sjrsten 
II, on the other hand, is less enplex in hardware and 
is an all digital system. The synchronisation between 
transmitter and receiver is achieved by using 5'MEz 
Rubidium clock from which all other desired frequencies 
are synthesised, 15 Bit shift register generates the 
PH sequence which is transmitted at 10 mega baud rate 
after phase modulation. The final carrier frequency is 
2100 MHz . 

At the receiver of Rake I, Pig. 3.4, the received 
signal is heterodyned with 2160 MHz to get IF of 60 MHz, 
which is then resolved into inishase and quadrature 
components by product demodulation* These are then 
correlated with the x shifted version of the reference 
PH sequence at 15 Rake taps having two correlators at 



:UB1DIUM 

CLOCK 


F16.3-3 TRANSMITTER 












GO MHz 



FIG. 3-4 RECEIVER FOR RAKE 
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each tap one for inphase and one for quadrature conponent, 
The correlation is achieved hy using current switch which 
is the heart of the system. In this system the signal 
processing is done throughout in the analog node. The 
expected dynamic range is 80 db. 

At the receiver of Bake II, Pig. 3 . 5 , after de- 
modulation the quadrature and inphase components of the 
signal are obtained in binary form. The components are 
then cross correlated with the t shifted version PN 
sequence by bit by bit nodulo -2 addition and counting 
at the 15 Bake taps, Barrow has shown that the counter 
output gives minimum error estimate of p (t, t) under 
low SBB conditions of the received signal. Hence G-auss- 
ian noise is added deliberately at the IP stage. The 
dynamic i*ange of this Bake System is about 20 db. 

The overall system characteristics are given 
as under s 


(a) 

Operating frequency 

- 2,1 GHz- 

(b) 

Transmission Technique 

- AM 

(c) 

IF Prequency 

- 60 MHz 

(d) 

Bounding signal 

- PH sequence of 



period (2^^ - 1) 



bits corresponding 



to 3.2767 nillisecs 
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IF NOISE 



FlG.3-5 RECEIVER FOR RAKE n 
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(e) Method of signal nodulation 

(f) Signal "bandwidth 

(g) Multipath resolution 

(h) Mode of propagation 

( i ) Link range 


-PS keying at 
10 M baikls 

- 10 MHz 

- 0,1 nicrosecs, 

- Overland scatter 

- 330 KMS 
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CHiiPTEIl 4 

SIMUMIION MODEIi OP THE THOPOSPHBRB 
4«1 Doppler Shift Tropospheric Scatter Model 

A doppler shift nod el suggested by Birkeneier [3] 
is considered in this work. Birkeneier has suggested 
that the scattering nechanisn can be described as consi- 
sting of a la.rge number of discrete scatterers confined 
to one or more relatively thin horizontal layers wherein 
the scatterers drift along with the wind. Such a discrete 
scatterer hypothesis lends itself to the digital com- 
putation of the total received scatter signal, at a given 
instant of tine, as the sun of the signal voltages recei- 
ved via each individual scatterer present in the antenna 
pa-ttem common volume, region. 

The doppler shift model shown in Pig. 4.1 locates 
the transmitter, receiver and scatterers in an X,Y,2 
co-ordinate system with the origin at the mid point of the 
chord line of length 2D that joins the transmitter and 
receiver. Assuming multiple scattering effects are 
negligible, then the ray path length L (X,Y,Z) is defined 
as the distance from the transmitter through the scatte- 
rer to the receiver. 

L(Z,T,Z) = [(D+Z)^ + + £(D-Z)^+T^+Z^I^ (4.1) 
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If a ^ earth radius is employed to accoun.t 
for refraction effects, then surfaces of constant ray 
path length appear as ellipsoids of revolution about the 
chord line (Z - axis) with the transmitter and receiver 
tas focii. A transmitted signal received via any single 
moving scatterer is shifted in frequency by a small amount, 
proportiona.! to the rate at which the scatterer crosses 
the constant ray path surfaces. The doppler shift f^ is 
given by 


-F _ 1 dl 

^ J dt 

Substituting (4.1) in <4. 2) 


(4.2) 


4 ul(l~4D^/L^) + VT + ¥Z 

" " (1-16 A L. 


(4.3) 


where 



along path wind speed 


Y 



¥ 



cross path wind speed 
vertical wind speed 


2 2 

Since, the factor (1 — 4D /L ) is small, the along path 
wind speed contributes very little to f^^. Therefore, 
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soatterers that nove in a direction perpendicular to the 
midpath plane and parallel to the constant ray path len- 
gth surfaces contribute very little to doppler shift 
of the received signal. The vertical winds are usually 
much smaller, than the horizontal winds. Essentially all 
received signal variation is due to cross path horizon- 
tal scatterer movement. Therefore, the layered, drift- 
ing with wind, discrete scatterer model of tropospheric 
scatter can be simplified for purpose of simulation 
through the assumption that scatterers exist only in the 
nidpath plane. 

Assuming that the crosspath wind speed is much 
greater than the vertical speeds, then (4,5) becomes 



2 

A 

2 

A 


TT ( ) 

V sin a 

V a 


(4.4) 


where a is the azimuthal angle of the scattering p>oint 
from the great circle plane. 

The propagation time delay, ^(2,7, 2) associated 
with each scatterer, relative to the shortest, or oarth 
grazing ray path is given by 
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%{Z,Y,7l) = I [1(Z,T,Z) - L (0,0,2^)] (4.5) 

where (0,0, Zg^) is the point where the grazing; ray 
intersects the Z axis as shown in Pig. 4.1 a-nd 0 is the 
speed of light. The surfaces of constant path length 
can also he described as surfaces of constant tine delay. 
Then, the doppler shift f^j can he written in terns of 
delay t, given hy 

-S- 

where f^ = frequency of the transnit ted signal. 

The connon volune can he divided into concentric 
shells of constant path delay increments as illustrated 
in Pig. 4.2. At nidpath plane for Z = 0, fron (4.5) 

1(0,T,2) - L(0,0,2^) 

T = — 

c 

2[(D^+T^+2^)'^ - (D^+Z(j^)^] 

0 

1-* 


or, 

or, 

or. 


(4.7) 
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Tlie equation (4.7) is a circle on T,Z plane § on the 
right hand side of the equation the only variable is 

i ' ■ ' ' 

the 3elay 'u* So by changing the values of constant delay 

T the radius of the circle on Y,2 .plane changes. For any 

'! ■ 

constant delay t , if the radius of the circle is R 
then E is given by 

= [ -^ + (D^ + (4.8) 

■ By sinulating the link parameters and the comon 
voluhe geometry of Kanpur - Rainital troposcatter link, 
the constant delay shells are illustrated in Fig. 4.2, 

Belay increments of 0,025 p.secs is given in the- , simula- 
tion program. 

Ant enna 

4 . 2 Antenna Gain Pattern and/off Set Effects 

The signal voltage observed at a scatterer result- 
ing from a signal transmitted from a narrow be.an-width 
antenna is a function of the antenna radiation pattern 
through the angular relationship between the,- signal -ray- 
path to the scatterer and the bore— sight line of maximum 
gain of the transmitting antenna. A similar angiiLar 
relationship holds between the signal ray path from the 
scatterer to the receiving antenna and the maximum gain 
boresight line for ' this_antenna^^ Then- the ^gain. p)attern 
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effect upon received signal voltage for a single sca- 
tter er is given by 

(4-*9) 

wliere 

= angle between transmitted signal ray path 
to scatterer and boreright direction of the 
transmitting antenna 

= angle between received signal ray path from 
scatterer and the boreright direction of 
the receiving antenna, 

Ror this simulation work, a G-aussian circular 
gain pattern for both transmitter and receiver antennas 
is assumed. Then antenna gain pattern and off set eff- 
ects is given by a function (G) as under as explained 
at Section 2*3. 

G = Exp [- (Y - (4.10) 

4.3 Angular Denendance Punction 

For isotropic turbulence, if a single scattering 
layer exists in the common volume, its edge on view 
would be super imposed on the ellipsoidal cross section 
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as shoxna in Fig» 4.2, The scattering cross section 
a(^) -as a fiaiction of scattering angled is of interest 
since it deterninos hew the reflectivity of scattering 
structure varies with position. 

The scattering cross section per unit volaine cr(0) 
is the power scattered towards the receiver per imit 
volune per unit solid angle per unit incident power 
density, Tliis is given by (Tatarski [7]) 

2 — ^ 

^(‘t) = 2% sin ■p vp (E - ) (4,3.1) 


where y- = angle leetween the electric field intensity/ 
of the incident plane wave and the dire- 
ction of the receiver 

d^CE - Ea>(.) = spectrun of refractivity fluctuations 


The vector 


(K - eIX 


) is defined as shown in the 


Figure below. 



/ind the vector K = K 




has a mgnitude given by 


T" 
,1 a 


2E sin 


& 


( 4 , 12 ) 
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For isotropic turbulence nodel, based on Tatarski’s 
theory (7)y the spectrum of refractive index fluctuations 
due to turbulence is given by 

(K) = 0.033 On^ (Z)“^ (4,13) 

2 

where On = refractive index structure function coefficient 


11 


for Tatarski's model which is now gene— 


P = 3 

rally accepted 

This leads to a scattering cross section. 


a (f?) = 2% sin^/^ x 0.033 Gn^(2E sin (4. 


For snail scattering angles, O can be written in terns 
of the paths* chord length 2D, the height of the sca- 
ttering point IS and the cross path displacement of the 
scattering point y 


^ 1 
^ = tan ^ [ 




. (Y^+Z^) 

] + tan ] 

(D-X) - (D+X) 


= [ 


Y^ + i 


D 


for X 4< D 


14) 
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ff (®) a sin”^ [ 


a 


2 2 
Z +Y 


D 




a [ 


-p/2 

— 

D ^ 

p/2 


Z +Y - 


(4.15) 


For Tatarski’s model, the scattering cross a (€>) is 
given ty 

11/6 

a (&) a [ — 5 — 2 — ^ 

Y +Z - 

In this simulation work p = 2, ll/3, 5 are considered 
in obtaining the delay, doppler power spectra. Birkene- 
ier [3] has designated this scattering angle dependence 
as the angular dependence function. When a layer of 
some thickness is considered, Z in (4.15) is taken as the 
average height of the scatterers called Z^, then a (B-) 
is given by 

cr (^) aC— ^5 p] (4.16) 

y^+Zq ■ 

When the turbulence is anisotropic with the 
horizontal correlation distance greater than the ver- 
tical correlation distance, then the angular dependence 


is written as 
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2 

D ^ R 

o (3) a[ 22 2^^. A>1 

o 

model 

For Tatarski's / (4.17)':3‘becones 

2 11 

D -g- 

O ( ) ® [ o' 6' 0 ] 

A Y +Z_ - 


( 4 . 17 ) 


4.4 Ambiigtiity PTmction 

Measurement of the scattering* function of a time 

frequency - spread tropo -channel req_uires an input 

signal having an ambignity function that is very narrow 

in both time and frequency co-ordinates. Otherwise the 

measured impulse response of the tropo-channel will 

be smeared due to the ambiguity function of the input 

signal. The signal used in the Bake System^ a carrier 

that is phase shift modulated by a long sequence^ 

is an ideal probe signal. By choosing, very small length 

of a single modulating pulse (t ) and very large repe- 
(To) ° 

tit ion period /of the PF sequence, we can have the ambignity 
function reasonably close to an impulse. Therefore, 
the measured scattering function which is the convolution 
of the true scattering function and ambiguity function 
can be taken as the true scattering function. In other 
words, smearing of the ambiguity function is ignored. 
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In this ’work, the parameters and are taken 
as 0,1 microsecs and 5.2767 millisecs respectively, thus 
making the ambiguity function reasonably close to an impulse. 

In this simulation work,;bhe ambiguit37- function ( Af) is conside- 
red unit;/ by choosing the above values of the parameters t;o and !• 

4.5 Maisnitude Weighting Ftinctlon M (I) 

In order to account for the possibility of layers at 
different heights having different scattering strengths 
a magnitude "weighting function M (^), whose value is diff- 
erent but constant for each layer is used. However, M (^) 
is considered unity when single scattering layer is simulated 
in this work. For multiple scattering layers only, we have 
considered different weightings for different layers, 

4.6 Delay Shell - Scatterin<s: Volume 

¥e have seen that the doppler freq^uency is proportional 
to the off - great circle displacement (y — displacement), A 
first order estimate of the shape of the delay, doppler spe- 
ctra can therefore be made from the relative portion of 
scattering .'.volume available as a function of doppler fre- 
quency (or y). In Fig, 4.2 by introducing an hypothetical ; 

scattering layer at a height of 5 HMS and of thickness 400 Mtrs, ■ 
we can visualise the scattering area subtended by the 
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scattering layer for any delay shell. For example let us con- 
sider the scattered power density at Rake tap Ho. 4 for re- 
lative multipath delay of -0.2 microsec. The scattering power 
seen hy the Rake tap Ro.4 is due to scattering from the 
shaded area i.e. the area subtended by the scattering layer 
and constant delay shells of 0.175 microsecs and 0.225 micro- 
secs as shown in Fig, 4.2. The resultant scattered power den- 
sity as function of doppler shift is shown in Fig, 4.2 based 
on the area within a certain doppler frequency range. 

The area subtended by the scattering layer and the 
constant delay shells is computed as under ; 

(a) For the region ABC in Fig, 4.5» the area is pro- 
portional to Z axis length subtended by the delay 
shell having radii and 

For radius R^^, given by 

Zl =f(R^^ - Y^) 

' ' Area = SH 2 - 

= SH2 - T®) 

(b) For the region BCDB, area is given by 
Area = SH 2 - SH^ 


( 4 . 20 ) 


( 4 . 21 ) 
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(c) In the region DEF, area is again proportional to 
the Z axis length subtended by the delay shell 
For radius * ^2 gi'ven by 

= Y(R2^ - Y^) 

Area = - Y^)- SH^ (4.22) 


Z 



4.7 Discussion 

The delay doppler power spectral density simulated 
in this work is a function of the following factors discussed 
in brief in the foregoing sections % - 

(a) Scattering cross section or angular dependence 
function a (©) , given by 
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(Id) Antenna gain pattern G-, given by 

& = Exp - [(T - Y )^/Y ^1 

c m 

(c) Area subtended by the scattering layer and the 
particular Rake tap delay shell 

(d) Ambignity function A^, which is considered unity, 

(e) Magnitude weighting fimction M (X), which is 
considered unity for a single scattering layer 
simulations. For multiple scattering layers simu- 
lation, different integer values are considered for 
each layer as weighting function. 

For a single scattering layer, the amplitude of the delay 
doppler power spectral density is given by 


Amplitude = a (9) x G x Area (4.23) 

For a multiple scattering layers, the amplitiide of the delay 
doppler power spectral density as given by 

Amplitude = a {_&) z G z- Area x M {V) (4.24) 

is used for simulation work. 
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GHAITEIl 5 

SCATTBllING FIOTCTIOES ~ DELAY, DOPPLER POWER SPECTRA 

5.1 Delay. Doiopler Locus for a Single Layer 

The delay doppler power spectral density for relatire 
niiltipath delays of 0,05 microsecs to'0#95 microsecs are 
illustrated in Figs, 5.2 through. 5.25. In these figures 
we notice the parabolic trace that a layer produces in a 
scattering function. The parabola results from the fact that 
the multipath delay in the signal received from a point 
scatterer varies as the sq^uare of the cross path position, 
while the doppler frequency varies linearly with cross path 
position. This is shown as follows : 

For a layer at height 2, we wish to find the delay t 

as we move along the layer in the cross path direction Y. 

Since the constant delay surfaces are circles, x will not vaiy 

2 2 

along the surface defined by the radius + Y ) 

From equation (4.7) 

t:(0,Y,Z) = § [(D^+Y^+Z^)^ - (D^+Z^.^)^] 


(5.1) 
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2 ^ 

t( 0,Y,Z) = ^ [D(l+ 5 ) - ]3(i+ 

^ 1 D- ■ 


= §[d{h4( ^)] 



( 5 . 2 ) 


4 VT 

SuTostituting the doppler shift relation (4.4) f-n = » 

where 1^ is the grazing ray path length which is approximately 
equal to 2D 


f 


D 


2 7T 
A D 


or 


^dAd 

2 7 


(5.5) 


SulDStituting equation ( 5 . 5 ) ia equation ( 5 . 2 ), we hare 

2 


fn^VX^- D 7^2 

t( 0 ,Y,Z) = -2 + 

4 V^-G G*D 


2 




C*D 


By sulostituting the Kanpur - Hainital troposcatter link para~ 
meters and the geometry, we have 

f ^ 

t(0,Y,Z) = 2,8 + — - *0655 (5.4) 

^ 49.5 

where t is in micro secs 

V is in metrs/sec 

Z is in kms f 


is in Hz 
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A cross path ■wind velocity (V) of lh»8 netra/sec is assumed 
for the simiilation work. 

Equation (5*4) gives the delay - doppler locus for a single 
layer. 

5.2 Simulation of Scattering E-unctions 

The scattering function is a plot of the power spe- 
ctral density as a function of doppler frequency shift for 
each value of propagation delay Scattering functions are 
simulated based on the simulation model develqped in Chapter 
4 to illustrate the effects of antenna azimuth off set, antenna 
elevation off set, types of layers and power indes p of the 
power law spectra. 

Digital computer simiilation programs are written in 
Fortran IV source language which are processed by IBM-1800 
digital computer and the scattering functions are plotted 
by IBM-1627 plotter. The listings of these programs are 
given at Appendix A. A series of simulation results in the 
form of scattering f'unctions are obtained in Figures 5,2 
through 5.25. 

To facilitate the discussion of scattering functions 
results, the figures are divided into seven sets based on 
antenna off set angles, number of layers and thickness 
of scattering layer. 
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FOTE ! - The scale factor of relative power density has 

loeen changed for each fi,gure in each set to dis- 
play the lower spectral levels more clearly. 

5.2.1 Scattering functions ; Set - 1 

The scattering function results of this set corres- 
pond to the simulation set up shown in Fig. 5.1 . The si- 
mulation program listings for Fig. 5.1 an.d for Pigs. 5.2 
through 5.4 are given at appendix A-1 and A— 2 respectively. 

The salient simulation parameters are given as under : 

a) Antenna off set (1 ) = 0 

h) Single scattering layer at a height (SHI) = 3 EMS 

c) Scattering layer thickness (SH 2 -SH^) = 400 Mtrs 

d) P = 2, y, 5 

e) Isotropic layer (A = 1) 

f) G-aussian circular antenna pattern of 3 EMS ra- 
dius at the midpath plane 

The resulting scattering functions are illustrated in Pigs. 

5,2 through 5,4. The following oloservation are made t 

a) The more apparent effect of increasing the value of 
p is the increase in power received. This is because 
as p increases the scattering cross section increases 
governed by equation (4.16), 



CONSTANT DELAY SHELLS 



FIG. 5-1 
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OELAY-DOPPLER power SPECTRA'. FOR 

NO ANTENNA OFFSET Yc= 0 Kms 


SH, = 3-0 Kms 
SHz^. 3 4 Kms 
P ^11/3 

Scale-.AMMO^ amt 



FIG. 5-2 




OELAY-DOPPLER POWER SPECTRA: FOR 

NO ANTENNA OFFSET 

Yc = 0 
P =2 

SH,=3 KMS 

SH 2-34 KMS 

~2 

Scale: AMT=10X AMT 



Doppler Shift (fp) in Hz 


Fin. R-'i 





DELAY-DOPPLER POWER SPECTRUM 
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b) Even though there are scatterers present in the 
assumed scattering layer in regions illuminated 
by antennas, where the relative delay is greater 
than 0,4 microsecs, the level of signals received 
at these greater delays is reduced considerably by 
t he angular dependance function (4.16), 

c) Each of the figures of this set tends to illustrate 
the parabolic relationship of doppler frequency 
and multipath delay described by equation (5.4). , 

5.2.2 Scattering functions - Set - 2 

This set corresponds to - the simulation set up shown in 
Fig. '4.2. The simulation program listings for Figs. 5.5 
through 5.7 are given at Appendix (A- 2). The important 
simulation parameters are given as under s 

a) Antenna azimuth off set (Y ) = 3 EMS 

b) Single scattering layer at a height (SH^) = 3 KM3 

c) Scattering layer thickness (SH 2 -SH^) = 400 Mtrs 

d) p = 2, 5 

e) Isotropic layer (A = l) 

f) Gaussian circular antenna gain pattern of 3 EMS 
radius at the midpath plane. 

The resulting scattering functions are illustrated 
in Figs, 5,5 through 5.7. And the following observations are 


made s 



DELAY-DOPPLER POWER SPECTRUM: 
AZIMUTH OFFSET BY 3 KMS 


^♦TENNA 


Yc = 3 KMS 

P 2 
SH| = 3 KMS 

SH2=3-4KMS 

Scale*. AMT 
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:;-^^:trAT-DOPPLER POWER SPECTRUM'.ANTENNA 
AZIMUTH OFFSET BY 3KM5 

Yc = 3 KMS 
SH^ = 3 KMS 
SH^ = 3'4KMS 
P == 11/3 
A = r -5 
ScalG:AMT = 10XAMT 



FIG. 5-6 


DELAY-DOPPLER POWER SPECTRUM: ANTENNA 
AZIMUTH OFFSET BY 3 KMS 


Yc = 3 KMS 
P = 5 
SHi= 3 KMS 
SH2= 3'4KMS 

Scale: AMT=r1Cf^AMT 


Doppler freq in Hz(fD ) ► 


FIG. 5-7 

i.j.T. r;^FUP 
OEHTR/u- I.ISRaHY 




a) Tlie scattered power density clianges by antenna azi- 
muth. off set to a considerable extent. 

b) As the antenna off set is increased, signals recei- 
ved from the higher delay regions associated with 
the layer becomes more significant and signals 
from the lower delay regions are suppressed. 

c) Parabolic, relationship between doppler frequency 
amd multipath delay is justified 

d}The power density increases as the value of p incre- 
ases 

e)fhe level of signals received for relative delays 
greater than 0,6 raicrosec is reduced due to the 
angular dependance function 

5.2,3 Scattering functions ; Set - 3 

This set corresponds to the simulation set up, 
shown in Pig. 5*8* The important simulation parameters are 
as under : 

a) a) Antenna azimuth off set (Y ) = 3 EMS 

o 

b) Antenna elevation off set = 1,8 EMS 

c) Single scattering layer at a height (SH^) =s 4,6 EM 

d) Scattering layer thichness (SH 2 - SH^) = 400 Mtrs 

e) p - 2, 5 . 

f) Isotropic layer (A = 1) 

g) Gaussian circular antenna gain pattern of 3 EMS 
radius at the midpath plane 
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DELAY-DOPPLER POWER SPECTRUM ; ANTENNA 
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DELAY DOPPLER SPECTRUM ; ANTENNA 
AZIMUTH OFFSET BY 3KMS AND 
ELEVATION OFFSET BY 1-8 KMS 


Yc-= 3-0 KMS 
P = 11/3 
SHi=4-6 KMS 
SH2=5-0 KMS 
Scale: AMT=io* amt 



1234 56789 


Doppler freq in Hz(fD) 


FIG. 5-10 


DELAY-DOPPLER POWER SPECTRUM: ANTENNA 
AZIMUT' OFFSET BY3KMS AND ELEVATION 
OFFSET BY T8KMS 


Yc = 3-0 KMS 
P = 5 

SHi = 4-6 KMS 

SH2=5-0 KMS 
-1 

Scale: amt=10*amt 
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The resulting scattering fimctions are illustrated in 
Figs, 5,9 through 5,11. The observations given in section 
5,2,2 are valid. In addition the follo'wing are observed : 

a) lor the same scale of Figs. 5.5 - 5.7 and 5.9 - 5,11 

. ■ ve notice a considerable reduction in the amplitudes of 
power spectrum for figs, at 5,9 - 5,11. This is 
because of elevation off set of the antennas, 

b) power received in the Bake taps from 1 to 6 is zero 
(for relative delay from .05 microsec to 0,3 microsec) 
due to the increase in layer height from where the 
scattering is assumed, 

c) The power spectral peaks are flat due to the constant 
area subtended by the scattering layer and the 
delay shells. 

5,2,4 Scattering functions ♦ Set - 4 

The simulation parameters are given as under : 

a) Antenna azimuth off set (Y ) = —3 2MS 

c 

b) Single scattering layer at a height (SH^) = 5 FMS 

c) Scattering layer thickness = 400 Meters 

a) p = 2 , U 5 

e) Isotropic layer (A = 1) 

f ) Gaussian circular antenna gain pattern of radius 
3 EMS at the mid path plane. 
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FIG. 5-12 
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DELAY-DOPPLER POWER SPECTRUM: ANTE NNA 
AZIMUTH OFFSET BY -3 KMS 


Yc = -3 KMS 
SH, _ 3-0 KMS 
SH^ = 3-4 KMS 
P - 11/3 

Scale;AMT =10^fAMT 



FIG. 5-13 


DELAY-DOPPLER 

POWER SPECTRUM- ANTENNA 

AZIMUTH OFFSET 

BY -3 KM S 


Yc = -3 KMS 


SH, = 3 0 KMS 

SHj= 3-4 KMS 


P - 5 


Scale:AMT=10XAMT 


&»1 


l!l 


% -A 

<? -os' 


-8 -7 -6 -5 -4 

Doppler freq in 


FIG. 5-14 





DELAY DOPPLER POWER SPECTRUM ANTENNA 
AZIMUTH OFFSET BY-3KMS AND ELEVATION 
OFFSET BY 1*8 KMS 

Yc=-3KMS 
P = 5 

SHi = 4 . 6 KMS 
SH2 = 5'0KMS 

Scale;AMT=io^-)f AMT 




.IS 


<? 


-7 -6 


Doppler freq in Hertz 


FIG. 5-1 7 





DELAY DOPPLER POWER SPECTRUM ANTENNA 
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FIG. 5-15 
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DE L AY-DOPPLER POWER SPEC TRUM: ANTE NNA 
AZIMUTH OFFSET BY -3KMS AND ELEVATION 
OFFSET BY 1-8 KMS 

Yc=- 3KMS 
SHi = 4-6 KMS 
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FIG.5-16 


66 



CONSTANT DELAY SHELLS 
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5.2,6 Scattering fuactions ; Set - 6 

The scattering function results of this set co- 
rrespond to multiple scattering layers in the common voltame 
as shown in Pig, 5.18. Each layer shown in Pig. 5.18 has 
different degree of anisotrophy (A) and magnitude weiging 
function The simulation program listings of Pigs. 

5,19 through 5.21 are given at Appendix . (-A-2). 

The iraportant simulation parameters are given as under : 


a) /ntenna azimuth off set (Y^) = 3 2MS 
h) Scattering layer heights, Anisotrophy (A) 
weighting function (M(1)) 


Layer Po. 

layer height 

Anisotropy 

M(Jt) 

1. 

1.8 

1 

11 

2 

2.2 

2 

10 

3 

2.6 

3 

9 

4 

3.0 

4 

8 

5 

3.4 

5 

7 

6 

3.8 

6 

6 

7 

4.2 

7 

5 

8 

4.6 

8 

4 

9 

5.0 

9 

3 

10 

5.4 

10 

2 

11 

5.8 

.11 

1 
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FIG. 5-20 
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FIG. 5-21 
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c) Thickness of each scattering layer = 400 Mtrs 

i) P = 2, 5 

e) Gaussian circular antenna gain pattern of 3 EMS 
radius at mid path plane. 

The resulting scattering functions are illustrated in Figs. 

5.19 through 5.21. The observations given in Section 5,2,2 
are valid. In addition the following are observed s 

a) The increase in power spectral density is due to 
the increase in scattering volume subtended by the 
time delay shells at the multiple scattering layers* 

b) The decrease in power spectral density for rela- 
tive delays greater than 0,4 microsecs is due to 
the increase of anisotropy with height and decrease 
in scattering strength (mr.gnitude weighting 
function M(-l)), 

5.2,7 Scatterini^; functions ; Set - 7 

The simulation runs in this set are made to observe 
the effect of (i) degree of anisotropy (ii) the layer thickness 
on the power spectral density of scattering functions for the 
single scattering layer set up given in Section 5.2.2, 

5.2,7 (A) Run 1 


The important simulation parameters are given as under 
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a) Antenna azimuth off set (T ) =3 ZMS 

0 

b) Single scattering layer at a height (SH^) s= 3 KISS 

c) Scattering layer thickness = 400 Metrs 

a) p = y 

g) Anisotropy (A) =10 

The resulting scattering functions are illustrated in Pig* 
5,22. The following observations are made : 

a) T.cro doppler power spikes are noticed at the 
relative multipath delays of 0,1 and 0,15 micro- 
socs 

b) For relative multipath delays greater than 0,2 
micro secs, the power spectral density is negligi- 
ble due to the anisotropic angular dependence fun' 
ction as given by 


. 2„2 


A^Y +Z' 


11 

] T 


c) The effect of anisotropy is distinct, when com- 
pared with Pig. 5.6 which is simulated for iso- 
tropic layer (A = l). 


5.2,7 (B) Run 

The same simulation parameters given as in Section 
5.2.7(A) are used for this Itan also with anisotropy (A) = 20 
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FIG. 5-22 
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" •iV-DOl i'LER POWER SPEC TRUM : ANTENNA 
IZ'P'UTH OFFSET BY 3 KMS 
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FIG. 5-23 
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Tho rceultinr acrittcrinc functions are illustrated in Fig. 

Tr.i >1j 3 ;rvations given in Section 5,2TlU) are valid 
fer thivi run aloe. 

.".7 (C) ‘\u i. ~ 3 

rho ir-i ortant sinulation parameters are given as under ! 

n) /uitcnna azimuth ioff set (Y ) = 5 MS 

0 

! ^ scattering layer at a height (SH^) = 3 KMS 

c) ^icattoring j.ayer thickness = 800 Metres 

1 ) i . 

^ 1 

'Hi'C rorulting rcr.ttcring functions are illustrated in Fig. 5.24 
I'.i’i given in Soction 5*2,2 are valid. In addi~ 

rlon th'; frllov/lng are observed s 

n' Who'i c 'ciparod with scattering functions given in 

Fig, 5.6, the scattering functions obtained for this 
Itun uro flat due to small variation in scattering 
■K'olv.r.’.'"> with cross path distance Y. 

5*2.7 (D) Hun - 4. 

Tho sai'C oiinulation parameters given as in Section 
5.2.7(0) arc used for this Run also with scattering layer 
■‘■hicknosn of 1200 Metros* The resulting scattering functions 
arc illustrated in Fig* 5.25* fUe observations given in 
Section 3*2.7 (C) are valid for this run also* 
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CHAPTER 6 

6.1 Refractive Index Structure Function Coefficient C^^^ 

The structure function has been used by Tatarski [7] 
to describe the state of atmospheric turbulence. For fluctua- 
tions in the refractive index the structure function can be 
written 


where 


\ (2:*) = 2 ?^^^ lo » r » (6.1) 

) , called the (g^2) 

structure function coefficient 

= a number close to 0.25 for sustained turbulence 
Lo = Outer scale of the isotropic turbulence in 
the inertial subrange (ranges from 10 to i; 

100 Mtrs) 

= inner scale of turbulence (order of millimeter? 



the average vertical gradient of refractivity 


Tatarski [7] has shown that when the radar’s effective 
wavelength matches the eddy scales falling in the inertial 
subrange, then, the radar cross section of the isotropic re- 
fractive eddies is given by equation (4.14) and for p = 


we have 
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CJ 


= 0.39 X 2 


-11/3 . 2 




(^/2) (6.3) 


vhcre /\ ~ -f^ho trans:r.itted signal wavelength. 


= tho scattering angle 

/rom (6.4), it It evident that for a troposcatter link having 

certain A and ^ , the power received is a measure of C 

laclionomlo n >neured ie a complicated functil of 

pVGrourn (P), tompGr>ature (T), mixing ratio (q), potential 

i'er,;pc.r'iturc ('•radient ( and potential mixing ratio gradient 

( la.) ^ 




Prom equation (6,2) 


can be evaluated 


’Irom radioaordo data taking the outer scale Lo to be the 
uo.j2:ial value o." 100 meters and a^ = 0,25. 


° ^ libra luat ion of 0 ^^ from fiake Data 

The power of the received signal, scattered at 
height Z, is given by the radar equation. 


P = 
r 


Vt. 

4nR^ 


A 


4nR‘ 


TT f 




whore = received power in Watts 

P.|. = Transmitted power in Watts 


G. ,G = Transmitting and receiving antenna gains 
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V = Scattering volume in meters^ 

H = Ilidpath distance in meters 
A = Vfav ©length in meters 
a(") “ Scattering cross section per unit volume in 
meters**^ as given in equation (6, 3) 

= Pilling fraction of the common volume, 

Tho equation (6.^0 depends on the assumption that isotropic 

O 

turbulence exists where is measured. If not we can expect 
a discrepancy between the power received and the predicted 
by {6. This discrepancy is accounted by the assumption of 

anisotropy in the refractive structure. 

For comparing the received signal intensity with the 
Gn^ parcin’oter as a function of height, it is essential to 
detorminn the power per unit scattering volume as a function 
of height. For scattering layers, the power received at sero 
dopplor shift observed in a particular Hake tape output is 
associ.-tea vrith the height of the tap's delay shell at the 
eroat circle plane. By suhstitutlng (6.3) in (6.4). 

^r ” h 4 4 nR^ 4 It 

A numerical integration of the poteer density over that 
portion of tho common volume that is seen hy each correlator 
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i3 performed as given in [4] to find a satisfactoiy conver 
Sion of received power to The power density of the 

receive i oif-nal due to single scattering from a point 



C4ii)^ 


>' 11/3 




( 6 . 6 ) 


v;hcre 1’^ ) = Transmitting antenna pattern at the print 

(X,Y,Z) 

= Hocoiving antenna pattern at the point (X,T,2) 

= Distance from transmitting antenna to the point 
at (X,Y,Z) 

Rj 3 = Diotanco from receiving antenna to the point at 

(X,Y,Z) 

2 

Since only the zero doppler power is used to estimate the 
evaluation of (6.6) is done for a point on (X,Z) plane for 
which Y = 0 

2 

r\ -6, = transmitting antenna elevation 

I*-! " D -(- X rl > 'Pl 

angle 

2 

^ _ -/hr. i /h r. = receiving antenna elevation 

'ii-2“D~X f 2 * P 2 

angle 
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= 3^ H- (D + X)^ 

Op p 

=-■ 7 ■: {D - X) 

tM 


.y 





+ D?1 


'Tho rol‘'’tivc poucr denssity from a scattering at a fixed 
height ic obtained by considering only the factor 


[ 




o ’ p 

H, " Ho 


-ll/3 

Sin ("^/S)] 


and int ■■/'.r.'vtiag this factor from -X^ to X^, where is the 
alon,"; ivh;h diatanco from the origin of the common volume 
coi-rci-j ondlng to the ith correlator. 




h 


I 

-h 


h‘'' 


Sin 


11 

“3 


(S/2) dx 


(6.7) 


?1 find the r.,aativG total received power for each correlator 
it ir= noceana:ry to integrate (6.7) over the vertical distance 
seen by each correlator. It is also necessaiy to include 
in this integration, the correlators gain function W( / 

IIoiico tho relative zero doppler power from a correlator 

io given by 


(0) = / 2 w (z) Pj,(z) az 

O' -i 


(6.S) 


(Zj-Zj^) = thloloiess of the chooseji Hake tap delay shell. 


■'hero 
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A cos^put©!* program is developed as per listings 
given in Appendix (A-3) for computing the relative sero 
doppler power received at each correlator. In this program, 
fcr'ms.'.’itting antenna gain pattern (Jl^) on (Z,Z) plane is 
an earned to be Gaussian circular given by 


= e* 


(Z - z,)' 


m 


(6.9) 


r.'hcrc Z = the height of the bore sight at the midpath 

plane from 'Iransnitter - Receiver Chord line 

p 2 

Z, =) Variance of the Gaussian pattern which is 2.6 KMS 


for Kanpur - Hainital tropo scatter link. 


¥o have identical antennas at transmitter and Receiver. So, 


2 i^) is given by 


- (Z - ZJ‘ 

= e 


( 6 . 10 ) 


m 


So, the equation (6.7) becomes 


Pi(Z) = 


^i ® 


X 


(^/2) 


Zm2 


dx 


-X. 


[Z^ + (D+X)^] [Z^ + (D-X) ] 


\.’here ^ 

D = 165 KMS 


For Kanpur — ITainital troposcatter 


link. 
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k triar. 'ular correlator gain function is assumed in this 
pro-^rnr-.. The conputed relative zero doppler power of the 
corrclntors is normalized hy the maximum value which was 
foinid to bo in the oocond correlator, set at the relative 
..'Uiltipith delay of 0.10 microsecs. The results are tabulated 
in Table !• 

TABLE 1 


Correl ator 


Relative Power in dbs 


0.05 

0.10 

0.15 

0.20 

0,24 

0.29 


-"2,68957 

0 

-0.45325 

-1.04542 

-2.48470 

-4.33455 


Tho tabic indicates that with a conmon volume uniformly filled 

v/ith scatterers tho observed relative zero doppler power from 

correlator one would be 2.68957 dB less than that from corre- 

2 

lator two to convert it to a relative measure of 0^ . 

By knowing the power received in each correlator 

2 

we can calculate the value oi • 
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For the existing troposcatter link between Kanpur - 
irainital (CW signal) the long term median of the received 
power is -79 dBm. From this we can have rough estimate of 
the* value from the equation (6.4) given by 

Py =s Watts (for -79 dBm median) 

^t - 

■t. = Gr = (for G^ = Gr = 45 db gain) 

V s 5# 3 X 10^ metrs^ 

-A * ^ metres 

R = 165 X 10^ metres 


Substituting the above values in equation (6,*^), we can find 
out the value of cr(Z) as 

a(Z) = 1.3750 X lO"®'® 


Prou equation (6.3) for latarskl’s model 

a(Z) = 0.39 X 0jj2 ^ 

. , 1.3750 X 10'^-^ ^ ^ 


-11/3 0 2 Sln"^^/^(®/2) = 1.3750 x 10 


11/3 ^ 2^1/3 


n 


.39 


metres”^'^^ 


= 1.47 X 10 
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Vr’lucf^ of in the range of 10”^^ to 10 around a height 
of 1 Kh a’^poar to be consistent over the Indian Subcontinent as 
jcr ’:hc ccudios carried out by Mazumdar [8], 
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CHAPIBE 7 

couonisiolir 


7 .1 Applications of the Simulated Scattering- Fmiotions 


(a) Fron the scattering Fiinctions obtained in Figs » 5.1 
through 5.25, it is possible to obtain a reasonable estimate 
of tho average cross path wind speed, layer altitudes and 
layer thicknesses. Equation (5.4) gives the delay - doppler 
locus for a layer given by 


2,8 


!d,!L 

y2 49.5 


- ,0655 


(7.1) 


where f^ = Average doppler shift for delay t 
V = ZAverage cross path wind velocity 
Z = Average scatter height 

Froi'i Fig- 5-6, for the relative multipath delays of 0,2 iJ-secs 
Mivi 0.55 (loccs, the average doppler shifts observed from the 
scattering functions are 1.5 Hz and 4.6 Hz respectively. Sub 
ntituting the values of v and fjj in (7.1), ve obtain the 
followin; ■ two equations, which we can solve for ? and Z- 

2 

0.2655 = 4^ 
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2 

and, 0.6155 = ^ + 

^ 49,5 

■ ' V = 12,2 Metres/sec 

and, Z = 3.3 KMS 

Tho sinulatod values of V and Z are 11.8 metres/sec and 
3,2 ffilS respectively. So we notice the estimated valuer 
arc very close to the simulated values. 

From Pig. 5.6, it is observed that only two correlators 
having relative multipath delays 0.1 and 0,15 ^secs have zero 
doijplcr power. From equation (7.1) at great circle plane 
(fjj = 0), we can obtain the heights Z 2 ^1 co 2 ?responding 
to the delays 0.15 md 0.1 hseos. By solTing for Sg md 
in equation (7.1) » we have 

2^ = 3.268 EMS 

= 2.863 23© 

Thoroforo, the scattering layer thiotaess is given ty 
(Z, - Zj_) which is calculated as 405 Metres. So, it is 
aeon that the layer thickness estlaated from the scattering 
rmictlona is vo^ close to the simulated value whioh is 

400 Metres. 
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(b) ./c can find the degree of anisotropy of the scattering 

lay era from scattering functions, by knowing the power recei- 
“/ed nt the Rake receiver correlators. From Figs. 5.22 and 
5*2,», it io noticed that the tracks display the increased 
angular dependence by having zero doppler spikes. The an- 
gular dependence for anisotropic model is given by equation 
(4.17) i.o. is, 


Cf(0) [ 




AY 


4 " 2 


5 -^ 


F/2 


Th'-ri-d'oru, by knowing the scattering cross section o’(5) for 
■..ah CO ■'^relator, we can obtain the value of A as a function 
c hc'lglit (Z), And it is observed by Birkemeier [4] that the 
■'■n hajtropy increases with height, 

(c) A airiulation study of the type carried out here enables 
on>> bo study the effects of antenna radiation pattern and 
pointing angle on the scattering function. Figs, (5.1) th- 
rough (5.25) illustrate the effect of changing antenna 
of foot G upon scattering functions. This technique offers 
ri ncans of estimating the performance of communication 
equipment under various known conditions of atmospheric 
nt.ructurc, 

(<1) Reaults of this simulation program can be employed 

in oonjmotlon with experimental Rake system data to 
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r-rovide increased knowledge of the changing meteorological 
conditions of the troposphere* Obtaining experimental data 
under changing noteorological conditions will be veiy diffi- 
cult because vo don’t know when such conditions will occur, 
ihoreforo, simulation studies of these Infrequent events will 
s -.r*' U3 the burd.cn of continuously collecting the data. 

7 S^l r ’:e3tions 

In this simulation program uniform scatterer distribu- 
tion is used* However, various scatterer distributions such 
■ c r/iaosion, poisson otc can be employed to determine which 
' ntribution most closely represent the physical situation. 
..cd the effects of various foims of scatterers can also be 
.•.oud- cd if equations that describe the reflection, refraction 
' l absorption characterstics of the scatterers can be 


.rritton. 
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APPEimiX (A-1) 


C 

G 

c 

V/ 

G 

0 

n 

G 

r% 

G 


fj 

G 


THIS PHOOa/HI SIMULA.TBS THE GOHSTMT DELAY SHELLS OE THE 
SXPSHI IBITTAL 

.OPOSC/iTTEiH LINK BETWEEN EMPUR AND NAINITAL. THE DELAY 
' ' u Uvi ARE 

GOGCENrHIC CI2CLES ON Y-Z PLANE WITH INCREMENTS OP 0.025 tiSB 
i’-iE .'.SSUriED GAUSSIAN CIRCULAR ANTENNA PATTERN AND AN 
iriPOTH ;riCAL 

SCATTL tIilG LAYER IS SUPERPOSED ON THE DELAY SHELLS. 

LINA PA I 

0— VELOCITY OP IN KMS/SEC 

D— HAU-’ THE DISTANCE BETWEEN TRANSMITTER AND RECEIVER IN EMS 
y,. AAD 2A— Y AND Z COORDINATES IN MS OP THE CENTRE OP THE 

nr v;ls op Ti-U3 

r::iA PAnrsw. these coordinates keep changing deptbiding 


a 
0 
a 
c 
c 
c 
0 

1 1 JOB 
// FOR KLINE 
C 

c 

n 


..RTF: 
vrOR T.!E 

..ZliPiTiL.L AND ELIVATION ANTENNA OPP SETS. 

TMii ;iADIUS OF thr CIRCULAR ANTENNA PATTERN AT THE 
:itKr..TH IN KMS 

SUl— HCIGIIT OF THE SCATTERING LAYER IN KMS PROM THE 
TitA:sklTTBR 

■ RECEIVER CHORD LUB 

.9512— 3111 +THICKNBSS OP THE SCATTERING LAYER IN KMS 
Ztr— ^Mi/.ZING RAY HEIGHT IN KMS 


'THIS SUBROUTINE IS FOR DRAWING THE DOTTED LINES AT THE 
liElGiiT ZG, SHI AND 
S 


JH2 'indicating GRAZING RAY HEIGHT AND SCATTERING LAYER 

suhroutinb dlinb (z,yino) 

Y = *0 

CALL DFPL0T(1,Y,2) 

1 Y=Y+YINC 

CALL FPL0T(2,Y,Z) 
y=Y+yiNC ^ 

Call fplot(i,y,z) 

2 RETUiUI 

// po;{*':miN 
NOirpROCESS PiiOGiiAM 
ONE JO id INTEGERS 
100 S( PLOTTER, CAxRD) 
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30 


31 

32 

33 

101 


102 


105 

35 


DATA C /3,*B5/ 

DATA D,ZG / 165., 1.8 / 

CALL 3CALP( ,6, .6, ,0, .0) 

PLOTTiSIi ROUTINES FOR DRAWING THE Y-Z GRID LINES - 
CALL PGRID(1,.0,.0,1.,13) 

CALL RGRID(0,.0,.0,1.,13) 

RilADC 5,100) SHI, SH2 

.iEAD(5,10O)YA,ZA,R 

PLOr nUG THE DOTTED LINES STARTS HERB 

CALL DLIiCE(ZG,.25) 

CALL DLINB{SH1,.125) 

CALL DLIHB(SH2,.125) 

PLOTTING OF THE CIRCULAR ANTBMA PATTERN STARTS HERB 
C3ANG=:(2A-ZG)/R 

.UUJL=ATM( SqfiT(l.-CSANG ^.*2)/0SANG) 
AIMC={3.14159-ANGL)/100. 

IFEN=~2 
DO 30 1=1,201 
Yr;Y*.+H SIiI{AN0L) 

Z-=ZA-H COSfANGL) 

CALL FPLOT(IPEN,Y,Z) 

IPSN=0 

A I i G L=j li) G L+A INC 

NOMOO UiNS WRITING WITH THE HELP OF THE PLOTTER STA.RTS HERE 

HT=(3iI2-SHl) *.36 

rF(HT-.15)31,31,32 

CLiLL FCiL;R(7.,SHl+(SH2~SHl) *.2,HT ^67, HT,.0) 

GOTO 33 

CALL PC!Un(7.,SHl+{SH2-SHl) .15, .0) 

OOtrTIiiUE 

V/iaTE(7,101) 

FOR:l;T( » SCATTERING LAYER’ ) 

CALL FCfLJi{7-,1.35,.10,.15,.0) 

WRITS (7, 102) 

PORI'UTC’ GRAZING RAY HEIGHT’) 

Y =«15 

D0"?5 I==l»14 

CALL FGIL\R( -.6 , Y , . 10 , .15 , . 0 ) 

WRrps(7,i05)J 

R0RMAT(12) 

Y=Y+1. 

X=~.3 

DO 40 1=1.14 

CALL FC!Ua(X,-.5,.10,.15,.0) 

’fRITB(7,105)J 
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40 


10 i; 


104 


n 

* 


106 


C 

ICC 


0 


10 


4 


Vj 


99 


// >C£Q 
ociiinD 

h? 

j) • 0 
0.0 


Z-^X+l. 

UH(7,0,-1,.1Q,.15,.0) 

•; iT;E(7, 103) 

0:'!i''«VT( 'Y * DISTANCE PROM GREAT CIRCIS IR KM’) 

':ALL •'C!UR(-.8,5.,.10,.15,1.5708) 

WkirE(7,104) 

i'OR;ivT(*Z - HEIGHT ABOVE SURFACE IE KILOMETRES') 

C.;LL FO;LiR(8.,12.,.12,.18,.0) 

W RITE (7, 106) 

?O.I:LU'('COuSMT DELAY SHELLS') 

■fOMO AUMS WRITIIIG ElIDS HERB 

T0yi7-~HEL;.TIVE DELAY CORRESPONDING TO THE FIRST DELAYSHEL. 
TOUIil-- aBlATIYE DELAY INCREMENTS FOR THE SUBSEQUENT 
DSLi\Y SHELLS 

TO’JFV— HEL/.TIVE DELAY CORRESPONDING TO THE LAST DELAY SHE. 
.t'iAD( 3 , 100 ) TOUIV ,TOUrN , TOUFT 
fo:nlu’(3P8.3) 
rOUFV=TOIJFV+TOUIN ■». 5 

Iirr TALI SAT ION OF THE RBLATiyS DELAY TON 
TOi>^TOUIV 

.1 ;?IALi3ATI0N OP THE Z COORDINATE OF THE DELAY SHELL 


Z ZG 

DZar:3aRT(D*D+ZG ZG) 

CT-{ .5 *C ’»rOU*l.B-6+D2G) ^2-D*D) 

PLOfTTNC; OF THE FIRST DELAY SHELL STARTS HERB 
CALL FPLOT (l.SQRT(CT-Z ^) ,Z) 

Z I C = ( SQIiT ( CT ) -Z G) /200 . 

XI'SII=2 

DO 15 1=1,200 


5=Zt-ZDTC , ^ s 

JALL PPL0T(IPBI'I,SQRT(CT-.Z*Z),Z) 

iR.OrTING OF THE FIRST DEL/i-Y SHELL ENDS HERE 

Oi’LATIVE DELAY FOR THE NEXT DELAY SHELL IS SET HERE 


roij=rou+TuuiN 

'..F ( rOU--TOUF?)lO ,10,99 


STOP 

3:iD 

I'LAIN 


3.4 

3.0 

.025 


5.0 

1.0 
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APPENDIX (A-2) 


C 


f: 


c 


rt 

u 


a 


c 

// 

// 


THIS PROGRAM SIMUIATBS THE DELAY-DOPPLBR POWER SPECTRA 
OP THE 

PROPOSED EX^EfilMENTAI TROPOSOATTBR LINK BETWEEN KANPUR AETD 
HAINTTAL 

USUJ'J RAKR SYSTEM -IT IS ASSUMED THAT THE SCATTERERS ARB 
u..'i:i’aRMLy 

DISTHIBUTSD IN A DISCRETE SCATTERING LAYERS -GAUSSIAN 
CIRCULAR ANTENNA 
liADtAriOH PATTIRN IS ASSUIIBD 
LINK ILvilAIffiTBRS 

SIU-— HEIGHT OF THE SCATTERING LAYER IN KMS FROM THE 
I’lLLR-SMITTER-RECBIVER CHORD LINE 
SiiR— 3H1+THICKNBSS OF THE SCATTERING LAYER IN KMS 
Z ;»-.G.l,LZIiJG RAY HEIGHT IN MS -n.T 

D— -HALF TiOS DISTANCE BETWEEN TRANSMITTER AND RECBIYER IN KMo 
2ZSRO«-Mi5M HEIGHT OF THE SCATTERERS IN KMS GIVEN BY 
( jiil+ jH2 ) *0 5 

yC-.~AllT.;miA’AXIHUTH OFF SET IN KMS AT THE 
ynSQ—ALTENNA PATTERN VARIANCE IN KMSQ ON, Y— AXIS. 

YM3L IS 2-6 KMSQ FOR THIS LINK 
?OtJ— -HEL^VTIVE MULTIPATH DELAY IN MICROSE^ 

tllKS TAPS ARE SET-IN THIS WORK THE RAKE TAPS SARE SET 

Tu'a'iD°!'0U2S°m¥ER®MB OPPBE KJMIPATH BEI4TB 

.„,n. ..ORBB MT.H 

JO 15 

, FOR CTK 

I0C3( CARD, PLOTTER) 

OHB lOHD INTBGBiRS 

^'‘’‘'‘■‘S5ToofsL,8H2,ZS,A,B.Z2EE0,YC,rHSQ,T0 

■3(:LlC-SaiB‘loR OTMTHAL DBSSITT 

BISFLACEMMT OF 

ffisfiDCsJSjSOff, SOU'.rolSF.TDISI 

laSlf MOTH OF I0¥ AXIS n TBM OF HO OF FDISP 

IHCIfflNTS „„„ dOPPIBR FRBQ AXIS 

D?mV--B 0PP1M™MQ DEAWnfO IHOIIHED GRID 

nFra^®DOPPLEE FRBQ IHOEIMBHIS IK HBifflZ 
S(i:?™N2.PFV,BFRIK 


100 


101 



o o o o o 
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102 FORMAT ( 214 , 2'F8 *1 ) 

CALL SCALP(SCLP,SCLP,LPRIV, .0) 

m)isp=m)isp/scLP 

7DISP=VLISP/SCLP 

msp=VDispm 

!rHDSP=:HLISP m 

B=DPRIV 

C DRAWIRl OP irclutbd g-rid lines starts here 
DO 10 1=1, N2 
call PPLOT(-2,B,*0) 

GALL PPL0T(-1,B+THDSP,T7DSP) 

C DPlAWING op inclined grid lines ends here 

10 B=B+DPRIN 

C INTIALISTTION OP HORIZONTAL AND VERTICAL DISPLACEMENTS 

C CORRESPONDING TO INTIAL TOU 

HDISP=0,0 
VDISP=0.0 

CALL PPLOT(l,DFRrV, .0) 

DZG=SQRT (D*D+ZG*ZG) 

C R1 AND R2— INNER AND OUTER RADII OP THE CONSTANT DELAY 
SHELL CO.RRBSPONDING 
0 TO TOUl AND T0U2 RESPECTIVELY 

E1=SQRT( (0.15*T0U1+DGZ) ** 2—D 
7 R2=SQRT( (0.15*T0U2+DZG) ** 2-D*D) 

T0U=T0U2-0.025 

C Y AND YT — ARE THE ACTUAL AND PLOTTER PREQ CO-ORDINATES 
Y=0.0 
YT=HDISP 

CALL FPL0T(1,YT,VDISP) 

C BRANCHING OUT POR DIPPSRENT CASES POR THE PURPOSE OP 
COMPUTING THE 

C AREA SUBTENDED BY THE CHOOSEN DELAY SHELL AND THE SCATTERIN 
LAYER 

IP (R2-SH1) 4,4,71 

71 IP(R2-Se2)l2,12,S2 
82 IP(R1-.SH1) 22,22,72 

72 IP(R1-SH2)23,23,73 

DPRSV — DOPPLER SHIFT STARTING VALUE CORRESPONDING TO THE 
INTERSECTION POINT OP THE RADIUS R1 AND THE 
SCATTERING LAYER HEIGHT SH2 

DPRSl-POINT ON THE DOPPLER SHIFT AXIS CORRESPONDING TO THE 
INTERSECTION OP THE RADIUS R1 AND THE SCATTERING LAYER 
liEIGHT SHI 

C DPRS2 — POINT ON THE DOPPLER SHIFT AXIS CORRESPONDING TO THI 

INTERSECTION 

C OP THE RADIUS E2 AND THE SCATTERING LAYER HEIGHT SH2 


» 
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G DFRilV-— DOPPIER SHIFT BHD VAIUB^CORRESPQIDIirG TO THE 

IHTBHSECTION OP 

G RADIUS R2 MD THE SCATTERIIG LAYER HEIGHT SHI 

73 DPRSY=SQET(ia^-SH2*SH2) 

G FOR THE PURPOSE OF GOMPUTATIOI OF THE AMPLITUDE OP THE 
POWER SPECTRAL 

C DENSITY DOPPLER SHIFT FREQ AXIS IS DEVILED IN TO STEPS 
0.02 HERTZ 
N=DFRSV/0.02 
Y=FLOAT(N) 0.02+0.02 
YT=Y+HDISP 

CALL FPL0T(2,YT,VDISP) 

23 DFRSI=SQ RT(Rim-SHl*SHl} 

22 DFEEV=SQRT (R2*R2~SH1 ^Hl) 

DFRS2=SQET ( R2 *R2-SH2 ^H2 ) 

G BRANCHING OUT FOR THE- DIFFERENT CASES FOR COMPUTATION OP 
THE COMON VOLUME 

C AREA CORRESPONDING TO BACH CO.RRELATOR TIME DELAY 

6 IF(Y-DFRS1)1,74,74 

74 IF(Y-DPRS2)2,75,75 

75 IF(Y-DFREV)76,76,4 

76 A.REA=SQRT ( R2 ^2-Y*Y ) -SHI 
GO TO 5 

1 IF(Y-DFRS2)77,77,3 

77 AREA=--SH2-SQRT(R1*R1-Y*Y) 

GOTO 5 

2 AREA=SH2-SH1 
GOTO 5 

3 AREA=SQRT ( R2 *R2-I C ) -SQRT ( R1 m* J*Y ) 

5 GALL AMP(AEEA,D,A,Y,2ZER0,YC,YMSQ,VDISP,YT,HDJSP) 
YT=Y+HDISP 
GOTO 6 

12 DPREV=SQRT(R2*R2-SH1*SH1) 

IF(R1-SH1)19,19,14 
19 IF(Y-DFREV)7a,78,4 

78 AREA=SQRT(R2 R2-Y^tY) ,SH1 

CALL AMP ( AREA , D , A , Y , 2 ZERO , YC , YMSQ , VDI SP , YT , HDISP ) 

YT=Y+HDISP 

GOTO 19 

14 DFRSl=:SQRT(RbRL-SHl SHl) 

21 IF( Y-DPRS1)15,89,89 
89 IF(Y-DF:EIEV)79,4,4 

7 9 AREA=SQRT ( R2*R2-Y *Y ) -SHI 

GOTO 18 
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1 5 AREA=S QRT ( R2 ^2-T ^ ) -SQET ( El *E1-Y ^ ) 

18 CALI AEP(AEEA,D,A,T,ZZEEO,YC,yMSQ,¥DISP,YT,HDISP) 
YT=Y+HDISP 
&0T0 21 ' 

C GOMPUTATIOE POE PEBSBNT TOU BEDS 
4 R1=E2 

CALI PPIOT(2,9.0+HDISP,VDISP) 

'IOU2=TOU2+O.05 
HDISP=KDISP+0.25 
VDISP=7DISP+0.25 
IP(T0U2-1.0)7,7,80 
80 CALL EXIT 
END 

SUBROUTINE POE G0I4PUTING THE AMPLITUDE OP THE POWER 
SPECTRAL DENSITY 
^ EOR GALO 

' SUBPOUTINE AI-CP(AREA,D,A,Y,ZZEEO,YC,YMSQ,VDISP,YT,HDISP) 

TUTC=EXP(-(YC-Y) ** 2/ (YHSQ) 

SCATC=(D*D/(A*A*Y*Y+ZZBR0 *ZZBR0) ) **1.8335535 
AMT=SOATG %RBA*ANTG C.0E-5+VDISP 
CALL PPLOT(0,YT,AMT) 

Y=Y+0.02 
EETUEIT 
END 

// XEQ GTE 
*OCBND 
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APPMDIX (A-5) 


C THIS PROGRAM SIMULATES THE RELATIVE ZERO DOPPIER 

G POWER OE CORRELATORS. THIS GIVES THE RELATIVE MEASURE 

G OF CIT SQURE 

// JOB 
// FOR CTR 

IOCS ( CARD , TYPEWRITER ) 

ORE WORD IHTEGERS 
HORPROGBSS PROGRAI4 

C THETA — SGATTERIHG ARGLE IR RAD PARTS , GAMA- 3 DB AtlTERRA 

C BEAM-WIDTH IR HAD PARS, 

C T0U1ARDT0U2 COxRRESPORD TO FIRST CORRELATOR TIME DELAY 

D SHELL IR MICRO SECS 

READ (5, 101) THETA, GAMA, T0U1,T0U2 

101 FORMAT (4F8. 4) 

C D-TRARSMITTER-REOErVER MIDPATH DISTARCE IR XMS,ZG — 

C GR/iZIRG RAY HEIGHT IR 

0 EMS,ZMSQ-VARIARCE OF ARTERRA PATTER!? ER miS 

READ(5,102)D,ZG,ZMSQ 

102 FORMAT ( 3 F 84 3) 

C THE FAOTOR( SIR (THETA/2) **- 11/3 IS COMMOR FOR ALL GOMPDTi 

0 ARD ARGF IS 

C THE SMLL ARGLE ARGLB APPROZIMATIOR OF THE ABOVE FACTOR' 

ARGP=(THETA^.5) **(-3.667) : 

C ZC-HEIGHT OF THE MIDPOINT OP THE COMMOR VOLUIIE FROM 

C TRA1?SMITTER-RBCEIVER 

C CHORD 

ZC=D ^HETA ■*€ . 5 
DZG=SQRT(D O+ZG ^G) 

G ZIARD Z2 ARE THE HEIGHTS WHERE THE DELAY SHELLS OP BACH 

C CORREL/iTOR 

C CUT THE Z-AZIS),ARD THUS(Z 2 -Zl). IS THE COMMOR VOLUME 

C THICKRESS OF THE . ' 

C CORRELATOR 

Z1=SQRT(0*15^0U1+DZG) w2-D*D) 

31 Z2=SQRT( (0.15 T 0 U 2 +DZG) *'^-D^) 

C ZIRO IS THE IRCRIMERTAL STEP OF THICKRESS (Z2-Z1) 

C CORRESPORDIRG TO 

C EACH CORRELATOR FOR EVALUATIRG THE OUTER INTEGRAL 

ZIRC=(22-Z1)*0.05 

C TOU-EELATIVE MULTIPATH DELAY OF THE CORRELATOR 

T0U=T0U2-0.025 

C Zr4BAR IS THE MEAN VALUE OF 211RD 22 USED FOR THE PURPO.' 

C OF COMPUTING 
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THE LBHGTH OH POSITIVE AXIS OF THE COMH VOIUIffi 
GOR-HESPOHDIHG TOEACH 

CORREMTOR DBHOTED BY ZE.IEIGTH EXTEHIDS FROM -IB TOZE 
ZIffiAH=(Zl+Z2) 0.5 

DIFFEOTT FORMULAE ARE APPLICABLE FOR COMPDTIHG IE 
FOR ZMEAH ABOVE 
MD BELOW ZO. 

IF ( ZffilAH-ZO ) 21 , 21 , 2 2 

21 IB=(2.0 ZMBAH)/(THBTA-GAMA)-D 
GOTO 25 

22 IE=D-(2.0 ZMEA1T)/(THETA+GAMA) 

XING IS TxHE IHG?3:EHTAL STEP OF GOMMOH VOLUME LBHGTH 
FOR EVALUATING 
THE IHHER INTEGRAL 
25 XnTG=ZE^.05 

IHTIALISATIOH OF THE OUTER INTEGRAL 
Z 1=0.0 

INTIALISATION OP 2 TAEBN TO BE THE MID POINT OP THE FIRSI 

INCREMENTAL STEP 

Z=Z1+0.5*ZINC 

COiyiPUTATION OF THE OUTER INTEGRAL STARTS ' 

DO 50 1=1,20 
S=I 

F=F1*P2,¥HERE , F1=F2 ARE THE GAUSSIAN CIRCULAR ANTENNA C 
PATTERITS ABOUT 2 AXIS 
P=EXP(-2.0 ^Z-ZC)^2/2MSQ) 

COTViPUTATION OF THE CORRELATOR GAIN FACTOR, W2,T'JHICH VARIES 
ZERO TO ONE FOR THE FIRST HALF OF THE CORRELATOR THICMES 
AND FROM 

ONE TO ZERO FOR THE LATTER HALF 
IF (1-10)24, 24, 25 

24 WZ=S ^.1-0,05 
GOTO 26 

25 ¥Z=2.0-S^0.1+0.05 
INTIALISATION OF TPIE INNER INTEGRAL 

26 XI=0.0 

INTIALISATION OF THE X VALUE, TAEEN TO BE THE MID POUT 
THE FIRST INCREMENTAL STEP 
X=— XE+O • 5 

COMPUTATION OF THE INNER INTEGRAL STARTS 
DO 40 J=l,40 

RISQ AND R2SQ A.RE THE SQUBRES OF THE DISTANCES OF THE 
SCATTERING POINT 

FilOM Tfl':'; RECEIVING ANTENNA AND TRANSMITTING MTTDNNA 
R-jSPECTIVBLY 


99 


RISQ=Z ^+(D+Z) **2 
R2SQ=2'^+(D-X) 

II=ZI+1.0/(RlSQ *S2SQ) 

40 x=x+ziNa 

C GOI'-IPUTATIOI OF THE HEER IFTEGEAL STARTS 

ZI=ZI+ZI ^Z*P 
30 Z=Z+ZI1T0 

C COT-IPUTATIOI OP THE OUTER IITTSGRAL BIDS 

0 ZEP— RELATIVE ZERO DOPPLER PO¥BR OP EACH GOIHIELATOR IE I 

ZLP=10 . 0 ^LOG( ZI1AEGP ) /2 , 3 
¥RITE(1,100)TOU,ZLP 
100 P0RffiAT(/10Z,P5.2,10Z,P8.5) 

C FOR COr'IPUTIEG ZDP OF THE EBIT CORRELATOR 

Z1=Z2 

T0U2=T0U2+0.05 

C AS THE VALUE OP Z2 OORRESPOEDIEG TO T0U2 GREATER THAE 

C 0.325 PALLS 

C OUTSIDE THE COMMOE VOtUME THE COJIPUTATIOE BEDS 

IP ( T0U2-0, 35) 31, 31, 52 
32 CALL EXIT 
BED 

//XEQ CTR 
CCEED 

0.0388 0.0166 0.025 0.075 

165.0 1.8 2.6 
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